function rest = randucoudw(t,~,T,a,b,m,sig0,alphaK,betaK,sig,s)


sigmavalue=[sig.sigma0(t,T,a,m,sig0),...
    sig.sigma1(t,T,a,b),sig.sigma2(t,T,a,b),sig.sigma3(t,T,a,b),sig.sigma4(t,T,a,b),sig.sigma5(t,T,a,b),...
    sig.sigma6(t,T,a,b),sig.sigma7(t,T,a,b),sig.sigma8(t,T,a,b),sig.sigma9(t,T,a,b),sig.sigma10(t,T,a,b),...
    sig.sigma11(t,T,a,b),sig.sigma12(t,T,a,b),sig.sigma13(t,T,a,b),sig.sigma14(t,T,a,b),sig.sigma15(t,T,a,b),...
    sig.sigma16(t,T,a,b),sig.sigma17(t,T,a,b),sig.sigma18(t,T,a,b),sig.sigma19(t,T,a,b),sig.sigma20(t,T,a,b),...
    sig.sigma21(t,T,a,b),sig.sigma22(t,T,a,b),sig.sigma23(t,T,a,b),sig.sigma24(t,T,a,b),sig.sigma25(t,T,a,b)];
mvalue=[s.m1(t,T),s.m2(t,T),s.m3(t,T),s.m4(t,T),s.m5(t,T),...
    s.m6(t,T),s.m7(t,T),s.m8(t,T),s.m9(t,T),s.m10(t,T),...
    s.m11(t,T),s.m12(t,T),s.m13(t,T),s.m14(t,T),s.m15(t,T),...
    s.m16(t,T),s.m17(t,T),s.m18(t,T),s.m19(t,T),s.m20(t,T),...
    s.m21(t,T),s.m22(t,T),s.m23(t,T),s.m24(t,T),s.m25(t,T),...
    s.m26(t,T),s.m27(t,T),s.m28(t,T),s.m29(t,T),s.m30(t,T),...
    s.m31(t,T),s.m32(t,T),s.m33(t,T),s.m34(t,T),s.m35(t,T),...
    s.m36(t,T),s.m37(t,T),s.m38(t,T),s.m39(t,T),s.m40(t,T),...
    s.m41(t,T),s.m42(t,T),s.m43(t,T),s.m44(t,T),s.m45(t,T),...
    s.m46(t,T),s.m47(t,T),s.m48(t,T),s.m49(t,T),s.m50(t,T)];
rest=zeros((alphaK+1)*(betaK+1),1);
rest(1:betaK+1)=0;
for ii=1:alphaK
    rest(ii*(betaK+1)+1:(ii+1)*(betaK+1))=sigmavalue(1:betaK+1)*mvalue(ii);
    
end
